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ELECTRON  COLLISION  SHIFT  OF  THE  LYMAN-ALPHA 
LINE  IN  H  AND  HE  + 

I.  Introduction 

>*The  interaction  of  plasma  microfields  with  radiating  atoms  and  ions 

causes  a  line  broadening  that  governs  both  the  the  half-width  and  shift  of 

the  spectral  line,  i.e.,  the  intensity  distribution  in  frequency  is  altered 

to  reflect  the  prese.  ce  of  the  plasma.  The  line  shifts  are  of  particular 

interest  both  for  their  intrinsic  value  in  terms  of  providing  insight  into 

the  fundamental  nature  of  the  underlying  physical  processes  affecting  the 

emission  of  radiation  as  well  as  its  diagnostic  value  in  the  study  of  high 

density  plasmas.  In  addition,  exact  wavelengths  of  spectral  lines  of  highly 

charged  ions  are  required  in  the  determination  of  line  opacities  and  level 

positions  for  x-ray  laser  transitions. 

Since  the  first  investigation  by  Berg,  et  al.^  on  the  effects  of  plasma 

interactions  on  the  frequency  shift  of  spectral  lines  it  still  remains  an 

interesting  and  challenging  problem.  This  is  due  primarily  to  its  N-body 

nature  for  which  only  approx  mate  solutions  can  be  found,  however,  in  low 

density  plasmas  the  mutual  interactions  of  perturbing  particles  may  be 

neglected  converting  an  intractable  problem  into  a  manageable  one,  making  it 

possible  to  evaluate  line  profiles  in  the  impact  approximation  via  the 

2 

Baranger  formalism.  This  theory  is  based  on  the  knowledge  of  the  scattering 

matrix  for  collisions  of  electrons  and  the  radiator  in  the  upper  as  well  as 

in  the  lower  state  of  the  line.  In  the  application  of  Baranger's  theory,  the 

initial  per turber-radiator  correlations  are  generally  neglected,  and  the 

theory  leads  to  negative  (i.e.  red)  frequency  shifts  of  Lyman-alpha  lines  of 

hydrogen-like  ions'*  ,  contrary  to  experimental  results  which  yield  blue 

9 

shifts.  It  has  been  suggested  that  another  effect  which  is  associated  with 

initial  per turber-radiator  correlations,  namely  plasma  polarization,  may 

cause  blue  shifts  of  hydrogen-like  lines.  Plasma  polarization  effects  on  He+ 

lines  have  been  calculated  by  several  authors^  with  the  assumption  of 

charge  neutrality  outside  the  orbit  of  the  bound  electron  with  the  result  of 

blue  line-shifts.  However,  a  procedure  similar  to  the  quantum-mechanical 
'  1 0 

approach  of  Volonte  ,  but  with  inclusion  of  electric  charges  in  the  outer 
region,  leads  again  to  red  shifts  of  Lyman  lines. 

In  the  present  investigation,  we  concentrate  our  efforts  on  a  few  special 
problems  in  J,e  theory  or  tiie  irne-shi ft ,  mainly  the  sensitivity  of 
calculated  shifts  to  various  effects  taken  into  account  in  the  evaluation  of 
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scattering  matrices.  For  this  reason  ve  have  chosen  the  Lyman-alpha  lines  of 

hydrogen  and  of  ionized  helium  as  two  typical  examples  where  the  effects  can 

be  demonstrated  and  differences  between  a  neutral  and  an  ionized  emitter  may 

be  seen.  Our  approach  is  fully  quantum-mechanical,  and  the  scattering 

matrices  are  obtained  in  a  distorted-wave  approximation  with  exchange.  This 

approximation,  although  rather  simple,  is  convenient  for  our  purposes  because 

it  is  very  well  suited  for  the  calculation  of  various  effects  -  addition  or 

deletion  of  different  terms  in  the  scattering  matrix  can  be  performed  easily. 

We  calculated  only  the  electron  collision  contributions  to  the  shift. 

Effects  of  plasma  ions  will  be  neglected.  Our  basic  assumption  is  that 

per turber-per turber  interactions  are  weak  and  can  be  ignored,  and  therefore 

the  applicability  of  the  results  is  limited  to  low-density  plasmas,  where  the 

Debye  shielding  length  is  much  larger  than  the  dimensions  of  the  radiator, 
18  -3 

i.e.  for  Ng  <10  cm  .  However,  the  effect  of  shielding  has  to  be  taken 
into  account  even  at  the  low  density  limit  for  perturbing  collisions 
involving  the  An=0  transitions  and  it  will  be  separately  discussed  in  par. 
II. 4. 

The  quantum-mechanical  results  are  then  compared  to  line-shifts  derived 
from  considerations  based  on  the  electrostatic  interaction  of  the  bound 
electron  with  the  plasma.  It  will  be  shown  that  at  the  limit  of  weak 
interaction,  both  results  are  identical,  but  that  neither  in  hydrogen  nor  in 
He+  has  this  limit  been  reached. 

II.  Calculation  of  the  Lyman-alpha  line  shift  using  the  scattering  matrix 
formalism. 

The  shift  of  the  Lyman-alpha  line  of  hydrogen  and  ionized  helium  produced 
by  electron  collisions  will  be  calculated  with  the  assumption  that  the  impact 
approximation  is  valid,  which  means  that  results  will  be  applicable  to  low- 
density  plasmas  only.  We  will  also  assume  that:  1)  the  per turber-per turber 
interactions  are  small  and  can  be  neglected,  2)  the  wave  functions  of  the 
bound  electron  are  not  changed  by  plasma  interactions,  and  3)  we  will  ignore 
the  effect  of  positive  ions.  No  uniform  positive  background  will  be  included 
in  the  calculations,  and  the  fine  structure  splitting  of  atomic  levels  will 
be  neglected.  All  calculated  shifts  will  be  normalized  to  the  electron 

_3 

density  N  =1  cm  .  Unless  indicated  otherwise,  all  quantities  throughout  the 
paper  will  be  given  in  atomic  units  (e=h/2n=m=l ) . 
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Within  the  validity  of  the  impact  approximation,  the  electron  collision 

profile  of  a  spectral  line  is  determined  by  electron  scattering  on  the  upper 

and  the  lower  level  of  the  line.  Let  us  denote  by  £  and  £  the  orbital 

quantum  number  of  the  bound  and  the  colliding  electron,  respectively,  and  by 

a  subscript  p  or  s  all  quantities  associated  with  the  scattering  on  the  upper 

or  lower  level.  The  scattering  on  the  upper  level  n^f^  is  characterized  by 

elements  of  the  scattering  matrix  SCa^,  a^),  where  the  scattering  channels 

a  ,  a  are  defined  by 
P  P 


a  h  n  £  k  £  L  S  , 
P  PPPPPP 


$  t t  r  f 

a  s  n  £  k  £  L  S 
P  PPPPPP 


t 

kp  (k^)  =  momentum  of  the  colliding  electron  before  (after)  the  collision, 

LpSp  =  total  orbital  and  spin  angular  momentum.  Similarly,  for  the 

scattering  on  the  lower  level  n  £  we  have 

s  s 


a  =  n  £  k  £  L  S  , 
s  s  s  s  s  s  s 


f  t  _  t  f 

a  =n£k£'LS 
s  s  s  s  s  s  s 


For  the  Lyman-alpha  line, 


n  £  s  2p,  L  =  £  ,  £  +1,  S  =  0,  1, 
P  P  P  P  P  P 


n  £  s  Is,  L  =  £  ,  S  =  0,  1. 
s  s  s  s  s  ’ 
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According  to  Baranger  ,  the  angular  frequency  shift  Aw 
transition  n^p  -*  ngs  is  proportional  to  the  electron  density 
expressed  in  terms  of  diagonal  elements  of  the  scattering 
and  S(ag,  ag)  by 


of  the  radiative 

N  and  it  can  be 
e 

matrix  S(a  ,  a  ) 
P  P 


3 


A<w  - 


12 


(2Lp  +  1) 


n  N  k-1  £(k)  E  t 

e»Lp»s 


c(2ST  +  1)  Im[l  -  s(ap,  ap)  S*(«s,  aj]  dk 


(1) 


with  k=k  =k  ,  S  =S  =S  ,  t=C.  =i  .  f(k)  is  the  momentum  distribution  function 

ps  ps  p  S 

of  plasma.  electrons  which  is  represented  by  a  Maxwellian  distribution  in  the 
present  calculation. 

T 

For  each  value  of  l  and  S  (and  consequently  for  each  matrix  element 

S(a  ,  a  )),  the  sum  over  L  has  three  terms  with  different  elements  S(a  , 
s  s  p  p 

ap).  A  schematic  example  of  diagonal  matrix  elements  S(ap,  ap)  for  total 
angular  momenta  Lp=l,2,3  and  one  particular  value  of  is  shown  on  Fig.  1. 
Elements  with  1=2,  corresponding  to  different  terms  in  the  sum  over  Lp,  are 
indicated  by  X. 

It  is  convenient  to  rewrite  (1)  in  terms  of  transmission  matrices  T. 
From  the  relation  T  =  1  -  S  it  follows 


1  -  s(v  ap)s*(V  as)  =  T(v  V  +  T*(v  <XS)  "  T(v  ap)T*(V  as)  (2) 


and  the  total  line-shift  is  thus  decomposed  into  two  parts:  contributions 

from  direct  terms  T(a  ,  a  )  +  T*(a  ,  a  )  and  from  the  interference  term 
*  p  p  s’  s' 

-T(a,  a  )  T  (a  ,  a), 
p  p  s7  s 

II. 1  Evaluation  of  the  T  matrix  for  ionized  helium  and  hydrogen. 

In  the  evaluation  of  diagonal  elements  of  the  transmission  matrix  T  it  is 
generally  necessary  to  take  into  account  coupling  of  several  scattering 
channels.  However,  in  the  case  of  the  Lyman-alpha  line,  the  energy 
difference  between  the  n=l  and  n=2  levels  in  hydrogen  and  He+  is  much  larger 
than  the  average  kinetic  energy  of  plasma  electrons  at  the  temperatures 
considered  in  our  calculation,  so  that  only  a  very  small  fraction  of 
electrons  scattered  from  the  Is  level  are  affected  by  coupling  to  other 
channels.  Therefore,  we  simplify  our  procedure  by  omitting  all  other  levels 
in  the  evaluation  of  the  matrix  elements  T(ag,  ag).  On  the  other  hand,  we 


have  included  levels  Is,  2s,  2p,  3s,  3p,  3d  in  the  T  matrix  for  the 

evaluation  of  T(a^,  a^),  but  we  have  ignored  all  elements  which  do  not 

involve  the  2p  level  and  consequently  have  only  a  small  indirect  effect  on 

T(a  ,  a  ).  The  structure  of  a  T  matrix  used  for  the  evaluation  of  diagonal 
P  P 

elements  corresponding  to  one  particular  value  of  and  S  is  shown  on  Fig. 
2.  Matrix  elements  taken  into  account  in  our  calculation  are  marked  by  X. 

All  matrix  elements  S(a,  a')  have  been  calculated  in  the  distorted  wave 
approximation  with  exchange  (DVX).  In  the  following  we  give  a  brief  outline 
of  the  procedure. 

12 

First,  we  calculate  the  p  matrix  with  elements  defined  by 


„<“••')  -  -2[(1-S„<xO  *„('  e  r  «'  lT)  Ro(' 

oo 

.v  ,  ff  I  r  t'  lt  )  R°(e  e  r  r) 


t  l  i'  i' 


E  *a(: 

X=1  M 


(-i)S  So  sx(?  i  r  e  J)  *ex(?  e  r  «')]  . 


e  and  V  are  the  orbital  angular  momenta  of  the  bound  electron  before  and  after  the 

T  T 

collision,  i  and  l’  correspond  to  the  colliding  electron,  L  and  S  are  the  total 

T  —  —  T 

orbital  and  spin  angular  momenta  of  the  system  with  L  =  L  or  L  ,  l  <L'  i'  L  ) 

—  —  T  P  s  i? 

and  g^((L  l  l'  t’  L  )  are  coefficients  given  by  Percival  and  Seaton. 

The  direct  and  the  exchange  radial  integrals  Rd  and  R^  are  defined  by 


4  -  ri)  Fke  (r2) 


(p  n'fjrJ  Fk,e,  (r2)  dtjdrj, 


re-  i  (  ' 

1  r< 

(  —  \  r 

\ 

Jpr;  riJ  Fke  lr2/ 

•!  X+l 

r 

[Pn'C';  r2)  FR,r  (, 

:lj  drl  dr2  • 

(5) 
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is  the  lesser  and  r^  the  greater  of  r^  and  r^,  P(nf;  r)  and  P(n'f';  r)  are  the 
radial  wave  functions  of  the  bound  electron  before  and  after  the  collision 
normalized  to  unity,  and  is  the  radial  function  of  the  colliding  electron 
before  the  collision  satisfying  the  equation 


[nf;  r)  =  P2jnf;  r^j  dr^  ,  r^  being  the  greater  of  r  and  r^. 


V(nf;  r)  represents  the  monopole  part  of  the  electrostatic  interaction  energy  of 
the  bound  a.-.d  the  free  electron.  Z  is  the  asymptotic  charge  for  the  colliding 
electron  and  Lagrange  multipliers  C,  C'  are  non-zero  only  if  e.=e.  or  i=i' , 
respectively.  They  are  determined  so  that  the  functions  P(nt;  r)  and  P(n'if';  r  ) 
are  orthogonal  to  F^  (r),  i.e. 


Jp  [nt ;  r)  F^  (r)  dr  =  Jpfn'r;  r)  Fke  (r)  dr  =  0  .  (8) 

The  asymptotic  form  of  F^  is 

Fke(r)  “  k  1/2  sin[kr  ~  ~J~  in  +  lt~  en(2kr)  +  argrfe  +  1  -  +  xj  .  (9) 

The  phase  shift  x  depends  on  k,  l,  and  on  the  level  nf,  but  it  is  independent  of 
T  T  a 

Radial  functions  F^,^,  of  the  scattered  electron  are  calculated  in  a  similar 

way . 


From  p  we  obtain  the  matrix  S  using 


1 


s'  =  (1  +  ifi)  (1  -  ig) 


(10) 


The  s' 


matrix  and  the  scattering  matrix  S  are  related  by 


S 


ix 

e 


(ID 


j  t  i  ^ 

e  being  diagonal  matrices  with  elements  e  a  .  The  diagonal  elements  of  the 
transmission  matrix  are  then  given  by 


2  i  t 

T(a,  a)  =  1  -  S'(a,  a)  e  a 


(12) 


In  the  present  version  of  the  distorted-wave  approximation,  the  channels 
included  in  the  p  matrix  are  uncoupled.  However,  transformation  (10)  introduces  a 
certain  amount  of  coupling  into  the  scattering  matrix  S,  and  consequently  the  line- 
shift  calculated  from  eq.  (1)  depends  on  the  number  of  levels  included  in  the  g  and 
S  matrices. 

In  the  summation  over  t  in  eq.  (1),  we  have  included  values  up  to  C  =  1 3 . 
Sufficient  convergence  was  achieved  for  the  case  of  He+,  for  hydrogen  we  have 
extrapolated  partial  contributions  to  the  shift  assuming  their  proportionality  to 
i  q  and  determined  q  from  values  for  (1=12  and  13.  The  extrapolated  contributions 
were  always  less  than  5%  of  the  total  shift. 

Values  of  T(a  ,  a  )  and  T(a  ,  a  )  have  been  obtained  for  10  different  energies 
\  2  P  P  s  s 

k  of  the  colliding  electron.  The  energy  mesh  was  determined  so  that  it  would 
give  the  best  possible  accuracy  of  the  integral  over  the  velocity  distribution,  and 
at  the  same  time  it  would  take  into  account  the  discontinuity  of  T(a  ,  a  )  at  the 

P  P 

excitation  threshold  of  the  n=3  levels.  No  resonance  effects  below  the  n=3 
threshold  were  taken  into  account  in  our  procedure. 

The  shifts  for  the  Lyman-alpha  line  have  been  calculated  for  electron 
temperatures  corresponding  to  2.0,  2.5,  3.0  and  4.0  eV  for  He+,  and  to  0.5,  1.0, 
1.5  and  2.0  eV  for  hydrogen. 
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If  the  diagonal  matrix  elements  of  S  are  much  larger  than  the  off-diagonal 
elements,  the  scattering  may  be  approximately  treated  as  a  one-channel  problem.  In 
this  case 


S(  a, 


a) 


2it| 

=  e  a 


(13) 


and  the  phase  shift  h  in  our  approximation  is  given  by 


t|  =  t  +  1/2  arctan(2p(a,  a)/{l 
a  a 


p2(a,  a)}] 


(14) 


according  to  (10)  and  (11).  Equation  (13)  indicates  the  importance  of  phase  shifts 

H  for  the  calculation  of  frequency  shifts  Ago. 
a 

The  accuracy  of  phase  shifts  h  obtained  from  (14)  may  be  illustrated  by  a 

a  14 

comparison  with  values  calculated  by  Oza  using  a  pseudo-state  close-coupling 
method  for  scattering  on  the  Is  state  of  the  He+  which  is  shown  in  Figs.  3-5. 
Oza's  method  is  substantially  more  accurate  than  our  simple  DUX  approximation,  but 
the  differences  of  the  two  results  represent  not  only  the  differences  of  the  two 
methods,  but  also  the  fact  that  his  phase  shifts  are  affected  by  coupling  to  other 
channels  so  that  (13)  is  generally  not  valid.  This  explains  why  the  smallest 
relative  differences  occur  for  £=0  (S-wave)  with  the  largest  phase  shift. 
Apparently,  the  major  contribution  to  the  phase  shift  comes  from  pure  elastic 
scattering,  therefore  (13)  is  approximately  satisfied  even  for  Oza's  values  of  the 
shift,  and  the  discrepancy  may  be  caused  just  by  the  difference  of  methods.  As  the 
elastic  scattering  decreases  with  increasing  t,  coupling  to  other  channels  becomes 
more  important  and  the  comparison  of  the  two  quantities  eventually  becomes 
meaningless . 

II. 2  Results. 


One  of  the  goals  of  this 
different  contributions  to  the 
series  of  calculations  in  which 
Aw.  A  list  of  all  calculations 


investigation  was  to  establish  the  importance  of 
total  line-shift.  Therefore  we  have  performed  a 
certain  terms  were  omitted  in  the  expression  for 
is  shown  on  Table  1,  and  results  of  the  distorted- 
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wave  approximation  are  presented  on  Figs.  6  and  7  for  He+  and  on  Figs.  10-12  for 
hydrogen.  In  both  cases  we  have  found  red  line  shifts. 

Curve  A  on  Figs.  6-13  represents  the  most  complete  DVX  calculation  with  all 
levels  ls-3d  included  in  the  S  matrix  as  described  in  par.  II. 1. 

For  He+  (Fig.  6),  the  omission  of  levels  3s,  3p,  3d  leads  to  curve  D,  which 

indicates  that  the  higher  levels  have  only  a  very  small  effect  on  the  frequency 
shift  of  Lyman-alpha  line.  If  only  Is  and  2p  levels  are  included,  the  results  is 
given  by  B,  and  if,  in  addition,  all  exchange  terms  are  omitted  and  only  elastic 
terms  retained,  the  shift  is  represented  by  C.  A  glance  at  Fig.  6  and  comparison 
of  A  and  C  shows  that  the  pure  elastic  scattering  is  by  far  the  most  important 
contribution  to  the  frequency  shift.  In  this  particular  case  Aw  depends  only  on 
phase  shifts  x  ,  because  g  =  0,  S'  =  1  and  consequently 


Aw  =  -nNe  k  Xf(k)  £  (2£  +  1)  sin2  (x  p&-  t  dk  , 


(15) 


where  xp^  and  x^  are  phase  shifts  associated  with  the  scattering  on  the  2p  and 
Is  levels,  respectively.  Furthermore, 


sin2 


(v  -  Ts<)  • 


sin2x 


p£ 


sin2x 


s£ 


Asinx  .sinx  „si 
p£  s£ 


in 


pe 


-  x 


s£ 


/ 


(16) 


The  first  two  terms  on  the  right-hand  side  of  (16)  correspond  to  the  direct  terms, 
and  the  third  term  to  the  interference  term  according  to  (2). 

The  importance  of  the  interference  term  for  He+  is  shown  on  Fig.  6  by  curves  J 
and  K.  J  was  obtained  in  the  same  way  as  A,  but  the  interference  term  was  not 
included.  K  similarly  corresponds  to  C.  It  can  be  seen  that  the  inclusion  of  upper 
levels  and  of  exchange  contributions  significantly  changes  the  magnitude  of  the 
interferenc  term  even  if  the  total  frequency  shift  is  affected  minimally. 

Fig.  7  demonstrates  the  effect  of  exchange.  Curves  L,  M,  N,  0  were  obtained  by 
ignoring  all  exchange  terms  in  the  g  matrix.  In  addition,  functions  F^  with  C=0,1 
were  not  or thogonalized  to  bound  wave  functions.  Curve  L  was  obtained  with  all 
levels  included.  A  comparison  with  A  shows  that  the  omission  of  exchange  decreases 
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the  total  shift.  Curve  N  represents  pure  elastic  scattering,  and  it  would  be 
identical  to  C  in  Fig.  6  except  for  one  thing:  functions  F^  used  to  calculate  C 
were  or thogonalized  to  bound  wave  functions,  while  those  for  N  were  not.  The 
resulting  difference  is  negligible.  Curve  M  corresponds  to  L  but  without  the 
interference  term,  and  0  similarly  corresponds  to  N. 

In  Fig.  8,  our  result  for  ionized  helium  (A)  is  compared  to  the  semi-classical 

3  4 
calculation  of  Griem  and  to  the  R-matrix  method  of  Yamamoto  and  Narumi. 

The  semi-classical  result  is  about  35%  lower  than  the  DWX  method  at  kT  =  2.5 
eV,  but  the  agreement  improves  with  higher  temperature.  However,  there  are  basic 
differences  between  the  semi-classical  procedure  and  our  method.  The  semi- 
classical  result  is  derived  solely  from  the  second-order  terms  in  the  S  matrix 
elements,  and  in  the  expansion  for  the  interaction  potential  only  dipole  and 
quadrupole  terms  have  been  taken  into  account.  In  the  DWX  method  the  main 
contribution  to  the  shift  comes  from  the  monopole  part  of  the  potential  and  higher 
order  contributions  to  the  S  matrix  are  due  to  the  unitarizing  relation  (10).  The 
£  matrix  has  only  first-order  terms.  Therefore  a  partial  agreement  between  the  two 
methods  appears  to  be  purely  fortuitous.  A  more  detailed  comparison  of  both 
methods  will  be  given  in  par.  II. 4. 

The  R-matrix  calculation  yields  results  which  are  almost  three  times  smaller 
than  the  DWX  approximation.  Yamamoto  and  Narumi  included  in  their  calculation  the 
same  number  of  levels  as  we  did  (ls-3d),  but  there  are  several  other  differences 
which  may  cause  the  discrepancy.  For  example,  they  use  the  same  basis  set  of 
functions  for  the  scattering  on  both  the  upper  and  the  lower  state.  According  to 
Yamamoto  ,  the  potential  for  the  scattering  on  He+  was  taken  equal  to  (in  our 
notation) 

r~l[l  +  exp(-21/3r)]  (17) 

as  compared  to  our  potential  2r-'*'  -  V(r)  with  V(r)  given  by  (7).  In  order  to 
evaluate  the  effect  of  such  a  potential,  we  have  repeated  the  calculations 
(including  all  levels  up  to  3d)  using  potential  (17)  to  generate  F^.  Thus  the 
differences  of  phase  shifts  t  -  disappear,  but  the  diagonal  elements  p(a,  a) 
were  supplemented  by  direct  monopole  terms  which  are  missing  in  the  DWX 
approximation  according  to  (3)  and  which  represent  different  elastic  scattering  on 
the  upper  and  lower  level.  The  result  is  shown  as  curve  S  on  fig.  9  and  it  proves 


that  the  form  of  the  potential  itself  cannot  explain  the  differences  between  our 

result  and  that  of  the  R-matrix  method.  We  suspect  that  the  evaluation  of  elastic 

contributions  in  the  R-matrix  calculation  is  responsible  for  the  discrepancy.  We 

have  performed  two  additional  calculations  to  substantiate  this  suspicion:  Curve  T 

was  obtained  in  the  same  way  as  S  (i.e.  using  potential  (17)),  but  without  monopole 

terms  in  the  p  matrix,  and  U  corresponds  to  our  original  calculations  (curve  A), 

but  with  all  phase  shifts  t  ^  set  equal  to  Tg(,  (the  result  is  then  independent  on 

the  actual  value  of  t  ).  In  both  cases  the  effects  of  pure  elastic  scattering  were 

a 

removed  from  the  calculation,  and  the  agreement  with  Yamamoto  and  Narumi  is  much 
better,  but  still  not  satisfactory  over  the  whole  temperature  range.  A  possible 

inadequacy  of  diagonal  R-matrix  elements  in  ref.  4  has  already  been  pointed  out  by 

„  .  3 

Gnem. 

5-8 

Observational  evidence  concerning  the  Lyman-alpha  shift  in  ionized  helium  is 
to  a  certain  degree  contradictory  due  to  experimental  difficulties.  In  contrast  to 
our  calculation  which  yields  negative  shifts,  experimentally  found  shifts  are 
positive  (i.e.  blue)  and  of  much  larger  magnitude.  At  the  present  time  it  is  not 
even  clear  if  this  discrepancy  is  caused  by  experimental  inaccuracy  or  if  it 
represents  incompleteness  of  the  theory. 

For  hydrogen  (Figs.  10-12),  the  sensitivity  of  frequency  shift  to  the  inclusion 
of  atomic  levels  in  the  S  matrix  is  much  larger  than  for  He+. 

The  smallest  shift  is  obtained  with  pure  elastic  scattering  without  exchange, 
i.e.  with  g  =  0  (curve  C,  Fig.  10).  The  shift  even  reverses  its  sign  and  becomes 
blue  at  low  temperature.  The  gap  between  A  and  C  is  much  larger  than  for  He+.  The 
largest  shift  is  the  result  of  inclusion  of  the  Is2s2p  levels  (curve  D).  Addition 
of  the  3s  and  3p  levels  (curve  E,  Fig.  11)  decreases  the  shift  and  so  does  the 
inclusion  of  3d  (curve  A).  If  4d  is  included  instead  of  3d,  the  shift  again 
becomes  larger  (G),  therefore  it  is  not  possible  to  estimate  simply  the  effect  of 
levels  with  n  >  3  by  extrapolation.  They  may  be  important  and,  in  contrast  to 

ionized  helium,  their  effect  may  not  be  negligible.  Calculations  without  the 

interference  term  are  represented  by  J  and  K  (Fig.  10).  J  corresponds  to  A  and  K 
to  C. 

A  comparison  of  results  without  exchange  is  displayed  on  Fig.  12.  As  in  He+, 

these  results  were  obtained  with  no  or thogonalization  of  F^  to  bound  wave 

functions.  The  omission  of  exchange  and  or thogonalizat ion  substantially  increases 
the  shift  (L.).  Curve  C  (Fig.  10)  and  N  (Fig.  12)  representing  results  of  pure 
elastic  scattering  differ  only  by  or thogonal izat ion  of  F^f  with  (=0  and  1  in  curve 
C.  For  hydrogen  this  difference  is  substantial  in  contrast  to  He+  and  it  is  caused 
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entirely  by  an  opposite  sign  of  the  partial  contribution  with  L^=l  in  the  DUX 
calculation  (curve  C).  The  effect  arising  by  the  omission  of  the  interference  term 
is  indicated  by  M  and  0  which  correspond  to  L  and  N,  respectively. 

The  semi-classical  calculation  of  Griem  is  compared  to  our  DUX  result  (curve 
A)  on  Fig.  13.  There  is  a  discrepancy  by  a  factor  of  five  between  the  two 
calculations,  and  the  disagreement  appears  to  be  due  to  the  differences  of  both 
methods.  A  discussion  of  the  discrepancy  and  a  comparison  with  the  experimental 
result  will  be  given  in  par.  II. 4. 

The  contributions  to  the  total  frequency  shift  Aw  of  the  Lyman-alpha  line  of 
He+  and  hydrogen  from  angular  momenta  L  are  shown  on  Figs.  14  and  15  for  two 
values  of  temperature.  In  He  ,  the  major  contribution  corresponds  to  L  =1,  for 
higher  angular  momenta  the  contributions  decrease  monotonously.  In  hydrogen,  the 
Lp=l  partial  contribution  is  positive,  while  all  others  are  negative.  The  opposite 
sign  for  ^=1  is  caused  predominantly  by  orthogonality  of  F^  and  P(nt),  and  it  is 
also  the  main  source  of  difference  between  line  A  and  L  in  Fig.  12. 

II. 3.  Shift  of  Lyman-alpha  derived  from  electrostatic  interaction. 


Several  authors  have  studied  the  Lyman-alpha  shifts  from  the  point  of  view  of 

las 
17 


electrostatic  interaction  of  the  radiator  and  plasma. ^  Status  of  the  theory 


prior  to  1978  has  been  reviewed  by  Volonte  .  The  basic  assumption  of  this 
approach  is  that  each  atomic  level  is  independently  shifted  by  plasma  interaction 
and  that  the  frequency  shift  of  the  spectral  line  is  equal  to  the  difference  of 
these  level  shifts.  The  level  shifts  are  therefore  derived  without  any  regard  to 
the  radiative  process  from  a  static  model  of  atom-plasma  interaction. 

It  is  easy  to  see  the  relationship  between  one  particular  form  of  "static  line- 
shift"  and  formula  (1).  If  all  values  are  small  compared  to  unity  and  only  pure 
elastic  scattering  significantly  contributes  to  the  line-shift,  Aw  is  given  by  (15) 
and  sin2(Tp^  -  t^)  by  (16)  without  the  interference  term.  Assuming  a  spherical 

charge  distribution  of  the  atomic  electron,  we  obtain  from  the  perturbation 

20 

expression  for  t fie  phase  shifts 


sin2(v  ■  xsJ  ~  -4JKv;  r)  -  v(nss;  r)]  Kc  (r)]  dr 


(18) 
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Potentials  V(nC;  r)  are  defined  by  (7)  and  F?  are  solutions  of  (6)  with  C=C'=0  and 


-1 


-i  1/2 

(i.e.  F^  =  k1  r  x  spherical  Bessel 


with  ?(Z  +  l)r  4  -  2V(nC;  r)  replaced  by  2Zr 
function  for  a  neutral  radiator  or  a  Coulomb  function  for  an  ion).  A  spherically 
symmetric  density  distribution  ng  of  mutually  non-interacting  free  electrons  moving 
in  the  field  of  a  positive  point  charge  Z  with  a  momentum  distribution  f(k)  can  be 
written  in  the  form 


ne  (r)  "  »Jk‘lf<k>'~2  E  (2*  *  l)  K  (r)]  dk 

with  lyn  ng  =  N  .  The  energy  W(nt)  of  electrostatic  interaction  of  a  bound 
electron  in  the  nt  state  with  free  electrons  described  by  (19)  is 

W^nfj  =  4ftjr^  ne  (r)  V^nt;  rj  dr  • 


(19) 


(20) 


The  expression  for  W(nt)  diverges,  but  the  difference  W(ne)  -  W(n'e')  has  a  finite 
value. 

Let  us  define  a  static  frequency  shift  Awstatjc  of  the  Lyman-alpha  line  as 


Aw  .  ..  =  W(2p)  -  W( Is)  . 

static 


(21) 


Then,  with  respect  to  (18),  the  static  shift  (21)  is  identical  to  (15).  It  should 
be  noted  that  for  neutral  radiators  with  Z=0,  formula  (19)  represents  a  uniform 
density  distribution. 

It  is  convenient  for  the  following  discussion  to  denote  by  Aw  .  _ .  T  the 

6  J  static  I 

frequency  shift  (21)  obtained  from  the  uniform  electron  distribution,  and  by 

„  the  shift  (21)  if  the  functions  F.  .  in  (19)  are  Coulomb  functions  with 
static  il  k£  ' 

Z*0. 

In  the  static  shifts  just  described,  the  effect  of  the  bound  electron  on  the 
free  electron  density  distribution  has  been  ignored.  This  effect  can  be  taken  into 
account  by  replacing  f£  in  (20)  by  F^ ,  which  are  solutions  of  (6)  (with  C=C'=0) 
and  which  are  identical  to  the  functions  used  in  our  distorted-wave  approximation 
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without  exchange.  This  procedure  leads  to  another  version  of  a  static  shift,  and 

the  result  will  be  referred  to  as  Aw  _  .  T,T. 

static  III 

In  a  real  physical  situation,  (18)  is  not  always  satisfied  and  the  interference 
term  is  not  negligible,  so  that  the  static  shift  may  be  substantially  different 
from  the  shift  calculated  according  to  (1).  Static  shifts  defined  by  (21)  are 
compared  with  the  DWX  result  (curve  A)  on  Figs.  9  and  13  for  He+  and  hydrogen, 
respectively. 

For  ionized  helium,  Aw  „  „ .  T  (curve  R)  is  much  smaller  than  /lu  ,  .  TT  (Q) 

static  I  static  II 

and  Aw  _  „ .  TTT  (P).  The  difference  is  obviously  caused  by  a  much  weaker 
static  III  '  j  j 

interaction  of  the  bound  electron  with  uniform  electron  background  that  with 

elecrLons  attracted  to  the  radiator  by  the  ion  field.  A  relatively  good  agreement 

of  P  and  A  is  probably  a  pure  coincidence.  On  the  other  hand,  the  agreement  of 

Aw  .  T  (curve  Q)  with  the  distorted-wave  calculation  without  exchange  and 
static  ll 

without  inelastic  terms  (curve  0)  is  rather  poor  in  spite  of  the  fact  that  both 
results  were  derived  from  the  same  expression  (15).  This  indicates  the  inadequacy 
of  (18)  for  the  replacement  of  phase  shifts  by  perturbation  expression. 

In  hydrogen  (Fig.  13),  A  is  much  smaller  than  Aw  .  „ .  ,  (R)  and  Aw  .  . .  ,TT 

(P).  The  phase  shifts  for  hydrogen  scattering  are  not  small  and  therefore  (18) 
cannot  be  valid. 

II. 4.  Comparison  of  the  semi-classical  calculation  and  the  DWX  method. 

3  16 

On  Figs.  8  and  13,  the  results  of  the  semi-classical  (SC)  calculations  ’ 
(circles)  for  the  Lyman-alpha  shift  are  compared  to  the  values  obtained  from  the 
present  DWX  method  (curve  A).  There  is  a  large  discrepancy  between  the  two  results 
for  hydrogen,  while  for  He+  the  agreement  is  much  better,  although  not 
satisfactory. 

In  the  semi-classical  calculation,  the  first-order  contributions  to  the  line 
shift  vanish,  and  the  shift  is  obtained  from  the  second-order  terms  in  the 
expansion  for  the  S  matrix  elements.  The  second-order  terms  do  not  contain  any 
monopole  contributions.  Moreover,  the  interference  term  in  the  expression  (2)  is 
ignored,  and  the  line  shift  is  equal  to  the  difference  of  level  shifts.  The  dipole 
contribution  to  the  shift  of  an  nt  level  from  electrons  with  velocity  v  and  impact 
parameter  p.  is  proportional  to 
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(22) 


(p.  v)  £_max(C',  l)  |  <n '  e  '  |  r  |  n£  >  |  ^  B  , 
1  n'C' 


where  n'C'  are  all  possible  final  states  produced  by  collisions  of  electrons  with 
the  radiator  in  the  state  nt,  and  B  is  the  shift  function  depending  on  p.  ,  v,  and 
E^-  -  En,-, .  Quadrupole  contributions,  which  are  much  less  important,  are 
proportional  to  a  similar  expression,  and  they  involve  matrix  elements 
<n' C '  | r^  |nC> • 

In  the  present  DUX  approximation,  the  matrix  elements  p(a,  a)  and  S(a,  a)  are 
calculated  from  the  first-order  perturbation  theory  and  contain  no  contributions  of 
type  (22).  In  order  to  obtain  second-order  contributions,  one  has  either  to  solve 
coupled  equations  for  the  scattering  problem  or  to  employ  the  second-order 
perturbation  theory.  In  the  following  we  use  a  simplified  form  of  the  latter 
procedure. 

Expression  (3)  for  the  elements  of  the  g  matrix  may  be  obtained  from  the  Kohn's 
variational  principle 


P(a,  a')  =  -2  £  f  Y*  (H  -  E)  Y  d? jdr,  ,  (23) 

spinsJ 

where  Y^  and  Y  ,  are  trial  wave  functions  of  the  total  system  (including  spin) 
corresponding  to  the  channel  a  and  a',  respectively.  Y^  has  the  form  (for 
simplicity  we  use  here  unsymmetrized  functions) 


Y  =  Q  r,1  F.  . 
a  a  2  k  t 
a  a 


w 


(24) 


with 


Q  =  E 


C  C  .T  ,  /0  .  /0  CI 
c  a  aLT  ^,1/2  1/2  S 


m  m 

C 

m  nT 


m  m  M, 


m  m 

s  s 


M 


$(c  mmlr^pfn  C  ;  r,)  m  (■■  ) ,  (25) 

(a  s)  1  (  a  a  1 ;  (  a  a  s) 
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1 

I 

S 


where  all  symbols  with  bars  refer  to  the  bound  electron,  $r^P  is  the  wave  function 
<if  the  hound  electron  including  spin,  and  *  is  the  angular  and  spin  part  of  the 
function  for  the  colliding  electron.  Improved  values  of  diagonal  elements  p(a,  a) 
may  be  obtained  if  the  trial  function  in  (23)  is  replaced  by  a  corrected 
function  Tcorr  which  also  contains  terms  representing  scattering,  so  that 


Crr  ■  vA  t  M  *2  V  *~i  Gk  ,e  ,M 

~  ~  or  ta  ~  ~ 


(26) 


a'  a' 


The  asymptotic  form  of  G  should  be 

Gk,e,  (r)  ~  p(“’  “')  k«'/2  C0S[ka'r  -  \  *«'"  +  Z*V  Cn  Kc' C) 


a'  a' 


+  ar*  r  K>  +  1  -  iZka')  +  Ta']  ' 

G  F 

In  our  simplified  procedure  we  replace  k  ,i  ,  by  p(a,  a')  k  , t  ,  which  has  the 
same  amplitude  at  infinity,  but  has  an  incorrect  phase.  Therefore  the  values  of 
matrix  elements  (23)  will  be  subject  to  errors  which  may  be  difficult  to  estimate. 
However,  with  this  replacement,  (23)  assumes  a  very  simple  form,  and  the  resulting 
expression  for  the  second-order  contributions  to  the  level  shift  contains  terms  of 
the  type  (22)  and  is  therefore  directly  comparable  to  the  SC  formula.  Assuming 
that  elements  p( a,  a')  are  small  and  retaining  only  terms  up  to  the  second  order  in 
p( a,  a'),  we  obtain  from  (23)  for  the  corrected  diagonal  matrix  element 


pCOrr(a,  a)  =  p( a,  a)  +  2  £  p2(a,  a')  .  (27) 

a'  *a 

l 

The  same  result  is  obtained  if  one  uses  antisymmetric  function  and  (27)  is 
therefore  generally  valid. 

t  On  Figs.  16  and  18  we  compare  our  previous  result  (curve  A)  with  line  shift 

obtained  from  matrix  elements  pCOrr(a,  a)  according  to  (27). 
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For  He  (Fig.  16),  the  result  of  inclusion  of  all  second-order  terms  for 

T 

angular  momenta  up  to  L  =13  is  represented  by  curve  B..  .  However,  the  sum  over  all 
t  1 

L  is  divergent  due  to  the  contributions  from  the  2p-2s  transition  with  AE=0.  If 

the  second-order  terms  due  to  this  transition  are  excluded,  the  shift  is  given  by 

curve  Aj .  The  effect  of  the  2p-2s  transition  cannot  be  properly  treated  without 

consideration  of  plasma  effects  on  electron  collisions,  because  contributions  from 

T 

distant  collisions  corresponding  to  large  angular  momenta  L  are  substantially 

21 

reduced  due  to  plasma  screening  effects.  In  a  recent  paper,  Griem  has  obtained 
values  of  contributions  to  the  Lyman-alpha  shift  from  An=0  transitions  for  He+  and 

17  _3 

H  for  the  electron  density  N  =10  cm  '  following  the  method  of  Boercker  and 
22  e 

Iglesias.  If  we  add  these  contributions  to  the  values  given  by  curve  (scaled 
to  N  =1  for  the  sake  of  comparison),  we  obtain  the  final  result  for  the  Lyman-alpha 
shift  represented  by  curve  C^. 

The  inclusion  of  second-order  terms  into  our  DVX  method  makes  the  disagreement 

3 

between  the  Griem's  original  SC  calculation  (Figs.  8  and  17)  and  our  present 
results  even  greater.  In  the  semi-classical  calculation,  the  main  contributions  to 
the  shift  originate  from  the  2p-3s  and  2p-3d  transitions,  but  the  total  effect  of 
the  second-order  inelastic  term  in  the  DWX  approximation  (difference  of  curves  A 
and  A^  on  Fig.  16)  is  much  smaller  than  the  total  SC  shift  (Fig.  17).  In  order  to 
find  out  if  the  discrepancy  is  just  a  result  of  differences  inherent  to  both 
methods,  we  have  performed  a  quantum-mechanical  calculation  which  follows  closely 
many  features  of  a  SC  procedure:  we  have  used  a  unitarized  Coulomb-Bethe 

approximation  with  inclusion  of  second-order  terms  according  to  (27),  but  we  have 
taken  into  account  only  contributions  from  the  2p-3s  and  2p-3d  transitions.  The 
result  is  represented  by  curve  on  Fig.  17.  If  our  results  are  extrapolated  to 
kT=5  eV,  the  agreement  with  the  SC  calculation  is  good.  The  disagreement  at  low 
temperature  may  be  caused  by  the  fact  that  in  the  SC  procedure  the  integration  over 
the  velocity  distribution  includes  also  non-vanishing  contributions  from  energies 
below  the  inelastic  threshold.  To  estimate  the  affect  of  such  contributions,  we 
repeated  the  previous  calculation  and  set  the  values  of  all  matrix  elements  below 
the  n=3  threshold  equal  to  their  corresponding  threshold  values,  thus  obtaining  an 
upper  limit  for  the  shift.  The  result  is  shown  as  curve  on  Fig.  17,  and  it 
demonstrates  the  drastic  effect  of  contributions  from  energies  below  the  threshold. 
From  the  results  for  He+  shown  on  Figs.  16  and  17  we  conclude  that  the  difference 
of  line  shifts  obtained  from  the  DWX  method  and  from  the  SC  calculation  is  caused 
by  several  effects:  elastic  monopole  contributions  which  represent  a  dominant  part 
of  the  DWX  shift  are  not  taken  into  account  in  the  SC  procedure.  On  the  other 
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hand,  the  inelastic  second-order  dipole  terms  are  much  larger  in  the  SC  method  than 
in  the  DUX  approximation  (agreement  with  the  Coulomb-Bethe  approximation  is  much 
better),  and  contributions  from  energies  below  the  inelastic  threshold  further 
increase  the  line  shift  in  the  SC  calculation,  while  these  contributions  are 
ignored  in  the  DUX  calculation.  Several  other  details  of  the  semi-classical 
calculation  may  also  contribute  to  the  discrepancy.  E.g.,  the  method  of 

uni tarization  (eq.  9  of  ref.  3)  makes  some  allowance  for  higher-order  perturbation 
theory  terms.  In  the  DU  method,  the  uni tarizat ion  is  achieved  by  a  different 
procedure  (eq.  10). 

Turning  to  hydrogen  (Fig.  18),  we  find  that  the  effect  of  second-order  terms  in 

the  DUX  approximation  is  greater  than  in  He+  (curves  B1  and  A1 ) .  Supplementing 
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results  shown  as  by  contributions  due  to  An=0  transitions  given  by  Griem  and 

scaled  to  N  =1,  we  obtain  curve  C.  which  is  in  a  good  agreement  with  the 
e  1  23 

experimental  value  of  Gruetzmacher  and  Uende  ,  from  which  the  ion  quadrupole 
contribution  according  to  Griern^  has  been  subtracted.  (The  result  of  ref.  23 

should  be  taken  with  caution,  however,  because  of  uncertainties  caused  by  the  J- 
splitting  and  by  the  fact  that  the  measured  width  is  greater  than  the  shift.)  As 
in  He+,  second-order  inelastic  contributions  in  the  DUX  approximation  (difference 
of  curves  A  and  A^  on  Fig.  18)  are  smaller  than  the  total  shift  from  the  SC 

calculation  (Fig.  13).  On  Fig.  19  we  again  compare  results  of  a  unitarized  Bethe 

approximation  (with  the  inclusion  of  second-order  terms  and  contributions  limited 
to  2p-3s  and  2p-3d  transitions)  and  semi-classical  shifts  from  An=l  transitions^ 
(circles).  In  the  SC  calculation,  the  contributions  to  the  shift  from  energies 
below  the  n=3  threshold  do  not  vanish.  These  contributions  represent  the  effects 
of  resonances  in  the  elastic  scattering  which  are  ignored  in  the  DUX  and  Bethe 
approximations.  If  the  contributions  from  energies  below  the  inelastic  threshold 
are  excluded  in  the  SC  calculation,  one  obtains  results  shown  as  triangles  in  Fig. 
19.  However,  this  result  still  corresponds  to  a  finite  value  of  excitation  cross 
section  at  threshold,  and  consequently  it  has  to  be  larger  than  the  result  of  the 
Bethe  approximation.  Figs.  18  and  19  for  hydrogen  therefore  lead  to  a  conclusion 
similar  to  that  reached  for  He+,  and  large  values  of  SC  shifts  in  comparison  to  the 
DUX  results  are  caused  mainly  by  larger  second-order  terms  and  by  non-vanishing 
contributions  from  energies  below  and  in  the  vicinity  of  the  n=3  threshold. 

III.  Summary  of  conclusions. 
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Frequency  shifts  of  the  Lyman-alpha  line  for  ionized  helium  and  hydrogen 

calculated  in  the  DWX  approximation  are  negative.  The  main  contributions  to  the 

line  shift  are  due  to  elastic  scattering  on  the  2p  level,  but  the  shift  is 

substantially  affected  by  the  inclusion  of  inelastic  terms  in  the  scattering 
matrix,  and  the  effect  is  larger  in  H  than  in  He+.  In  hydrogen,  contributions  from 
levels  higher  than  3d  (not  included  in  the  present  calculation)  may  not  be 

negligible.  Second-order  terms  in  the  scattering  matrix,  representing  coupling  of 
different  channels,  are  very  important. 

The  interference  term  increases  the  absolute  value  of  the  line  shift,  and  the 
effect  is  larger  for  H  than  for  He+. 

Major  contributions  to  the  shift  come  from  the  lowest  angular  momenta  of  the 
colliding  electrons. 

The  static  shift  of  the  line,  defined  as  a  difference  of  electrostatic 

interactions  of  the  bound  electron  with  plasma  electrons  moving  in  the  potential  of 

a  point  charge  Z,  is  larger  than  the  shift  from  the  DVX  method,  but  the  agreement 

for  He+  is  much  better  than  for  hydrogen. 

The  present  DWX  results  for  the  shift  of  ionized  helium  are  almost  three  times 

4 

larger  than  results  of  the  R-matrix  method  by  Yamamoto  and  Narumi  .  A  possible 

explanation  is  the  difference  of  elastic  terms  in  both  methods. 

The  difference  of  the  DWX  results  and  the  semi-classical  calculations  of 
3  16 

Griem  ’  is  attributable  to  several  effects:  elastic  monopole  contributions  that 
represent  a  major  part  of  the  DWX  shift  are  not  included  in  the  SC  calculations; 
inelastic  second-order  terms  are  much  larger  in  the  SC  method  than  in  the  DWX 
approximation;  and  inelastic  contributions  from  energies  below  the  n=3  threshold 
representing  resonance  effects  in  elastic  scattering,  were  not  taken  into  account 
in  the  DWX  method. 

+  5-8 

Experimentally  found  positive  shifts  of  the  Lyman-alpha  line  of  He  are  in 

sharp  contrast  to  our  calculation  which  yields  negative  shifts  of  much  smaller 

magnitude.  Observed  electron  collision  shift  of  the  Lyman-alpha  line  of  hydrogen 

is  in  agreement  with  the  DWX  result,  if  second-order  terms  are  included  and 

21 

contributions  due  to  the  2p-2s  transition  according  to  Griem  are  taken  into 
account . 

With  respect  to  previous  discussions,  the  most  reliable  result  of  the  present 

21 

investigation  (combined  with  partial  result  of  Griem  for  contributions  due  to  the 

2.p-2s  transition)  is  represented  by  curves  C1  on  Figs-  16  and  18.  It  applies. 

17-3 

however,  only  to  electron  density  Ng=10  cm  for  which  the  results  of  ref.  21 
have  been  obtained.  Further  improvement  in  future  calculations  should  be 
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accomplished  by  properly  taking  into  account  second-order  terms  in  the  scattering 
matrix  either  by  using  close-coupling  results  or  by  rigorously  applying  the  second- 
order  perturbation  theory,  and  by  taking  into  account  the  effect  of  resonances 
oelow  inelastic  thresholds. 
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Table  1 


Resti  Its 

for  the  Lyman-alpha  shifts  ; 

shown  on 

Figs.  6-13. 

(DUX  = 

distorted-wave  approximation 

with  exchange 

DU  = 

distorted-wave  approximation 

without 

exchange) 

Curve 

Levels  included  in  the 

Remark 

scattering  matrix 

A 

Is2s2p3s3p3d 

DUX 

B 

ls2p 

DUX 

C 

ls2p 

DU,  only  elastic  terms  included, 

F^  orthogonalized  to  P(n6)  with  6=6 

D 

Is2s2p 

DUX 

E 

Is2s2p3s3p 

DUX 

G 

Is2s2p3s3p4d 

DUX 

J 

Is2s2p3s3p3d 

DUX,  same  as  A  without  interference  term 

K 

ls2p 

DU,  same  as  C  without  interference  term 

L 

Is2s2p3s3p3d 

DU 

M 

Is2s2p3s3p3d 

DU,  same  as  L  without  interference  term 

N 

ls2p 

DU,  only  elastic  terms  included, 

F^  not  orthogonalized  to  P(n6) 

0 

ls2p 

DU,  same  as  N  without  interference  term 

S 

Is2s2p3s3p3d 

DUX,  F,  generated  using  Yamamoto  and 

4  15 

Narumi's  ’  potential 

T 

Is2s2p3s3p3d 

DUX,  same  as  S,  without  direct  monopole 

U 

Is2s2p3s3p3d 

DUX,  same  as  A,  with  all  x  „=x  „ 

'  p6  s6 

P 

static  shift  Aw  .  .  .  TTT 
static  III 

(from  distorted  waves  ) 

0 

static  shift  t  TT 

static  II 

(from  Coulomb  functions  F^ 

R 

- 

static  shift  Aw  .  .  .  T 
static  I 

(from  uniform  electron  background) 


terms 


Diagonal  matrix  elements  S(a^,  a^)  with  L^=l,2,3.  Elements  corresponding 

to  t =?  are  marked  hy  X. 


AlldVd  (aao)  N3A3 


mm 
mm 
mm 
mm 
mm 
mm 

mmm^mmmm 


AlldVd 

(N3A3)  aao 


■ 
m 

H 


T3 

C\J  CO  CO 


T— 

CM 

1 

+ 

Q. 

c 

Q. 

■D 

CO 

CO 

Transmission  matrix  T  (for  one  set  of  values  L 
calculation  of  diagonal  elements  T(ap»  a  ).  Element 
present  calculation  are  marked  by  X. 


ENERGY  (RYD) 

S-phase  shift  (radians)  for  scattering  on  the  Is  state  of  He  . 

14 

Solid  curve:  pseudo-state  close-coupling  calculation  by  Oza  ; 
dashed  curve:  present  result  from  the  DWX  approximation,  eq.  (14) 


ENERGY  (RYD) 

P-phase  shift  (radians)  for  scattering  on  the  Is  state 
Solid  curve:  Oza^S  dashed  curve:  eq.  (14). 


PHASE  SHIFT 


Au)(10' 


Lyman-alpha  shift  of  He  .  Symbols  are  explained  in  Table  1 

3 

Cirles:  semi-classical  calculation  of  Griem  ; 

4 

dashed  line:  R-matrix  method  of  Yamamoto  and  Narumi  . 
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Fig.  14.  Partial  contributions  to  the  Lyman-alpha  shift  of  He+  from  total  angulai 
momenta  L^.  Contributions  correspond  to  line  A  "Fig.  6). 

Solid  line:  kT=2  eV;  dashed  line:  kT=4  eV. 
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Aa>(10  'sec 


kT(eV) 


Fig.  16.  Lyman-alpha  shift  of  He  . 

A  -  DWX  calculation  (see  Table  1); 

Aj  -  same  as  A  with  inclusion  of  second-order  terms  according  to  (27) 
except  for  contributions  corresponding  to  the  2p-2s  transition; 

B.,  -  same  as  A  with  inclusion  of  all  second-order  terms  (contributions 

1  Y 

corresponding  to  2p-2s  transition  included  only  up  to  L  =13); 

21 

-  Aj  plus  contributions  from  2p-2s  transition  taken  from  Griem 
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-Aw  (10  'sec 


2  3  4  5 

kT(eV) 

Fig.  17.  Lvman-aipha  shift  of  He+. 

Dj  -  unitarized  Coulomb-Bethe  approximation  with  second-order  terms 
according  to  (27).  Only  contributions  corresponding  to  2p-3s  and  2p- 
3d  transitions  included; 

-  same  as  D^,  but  elements  of  the  g  matrix  below  the  n=3  threshold  set 

equal  to  the  value  of  corresponding  elements  at  the  threshold; 

3 

circles:  semi-classical  approximation  of  Griem  . 


39 


-Aa)(10  'sec 


Ao)(10~  'sec 
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kT(eV) 

Fig.  19.  Lyman-alpha  shift  of  hydrogen. 

-  unitarized  Bethe  approximation  with  second-order  terms  according  to 
(27).  Only  contributions  corresponding  to  2p-3s  and  2p-3d 
transitions  included; 

circles:  SC  approximation  of  Griem^  (An=l  contributions  only); 
triangles:  same  as  circles  without  contributions  from  energies  below  n=3 
threshold . 
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